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Elasto-incrtial turbulence is a new state of turbulence found in flows with polymer additives . The dynamics 
of turbulence generated and controlled by such additives is investigated from the perspective of the coupling 
between polymer dynamics and flow structures. Direct numerical simulations of channel flow with Reynolds 
numbers ranging from 1000 to 6000 (based on the bulk and the channel height) are used to study the formation 
and dynamics of elastic instabilities and their effects on the flow. The resulting mechanism of interactions 
between polymer dynamics and the flow helps resolve a long-standing controversy in the understanding of 
polymer drag reduction and explains the phenomenon of early turbulence, or onset of turbulence at lower 
Reynolds numbers than for Newtonian flows, previously observed in polymeric flows. Polymers also point out 
an interesting analogy with the forward and backward energy cascade in two-dimensional turbulence. 



FOREWORD 

This article is dedicated to our mentor, colleague and friend, Prof. Parviz Moin, in honour of his 60th birthday, 
and without whom this research would have never been undertaken. First, because our interest in polymeric flows 
started through Parviz' impulsion and guidance more than ten years ago, when research on turbulent drag reduction 
by polymer additives was initiated at the Center for Turbulence Research (CTR) in Stanford, and which project we 
had the privilege to be part of. Secondly, because years later the umbrella of the CTR and its Summer Program 
gave us the chance to revisit and further deepen this topic. The collaborative and incredibly fertile atmosphere of 
the CTR, which Parviz created and nurtured, was indeed a major contributor to the results presented here. Finally, 
and foremost, because Parviz' vision and pioneering work in numerical simulations of turbulence paved the way for 
so many future research in this field. This is very well exemplified by the results presented here, which demonstrate 
not only the capacity of high-fidelity computations to accurately simulate experimental observations, but also to offer 
an unparalleled insight into the physics of such complex flows. We are forever in debt to Parviz for his tremendous 
help, support and inspiration. 

I. INTRODUCTION 

Polymer additives are known for producing upward of 80% of drag reduction in turbulent wall-bounded flows through 
a strong alteration and reduction of the turbulent activity^. The changes in flow dynamics induced by polymers do not 
lead to flow relaminarization but, at most, to a universal asymptotic state called maximum drag reduction (MDR). 
The early quantitative description of MDR 2 was semi-empirical, yielding a correlation for the friction factor and the 
Virk log- law mean velocity profile. Recently, Procaccia et a/P proposed a theory that derives, in the limit of infinite 
Reynolds number, an asymptotic log-law for the mean velocity profile remarkably close to Virk's. However, at low 
and moderate Reynolds numbers, the existence of a logarithmic region in the mean velocity profile of MDR is not 
verified by currently available experiments and simulations 4 . 

The theory of MDR derived by Procaccia et al. describes the action of polymers as an eddy viscosity linearly 
increasing with increasing distance to the wall. Most noticeably, the theory assumes that "In contradistinction to 
the picture offered by de Gennes, simulations (performed by the authors) indicate that the energy never goes back 
from the polymers to the flow; the only thing that polymers can do is to increase the dissipation" (From Procaccia 
et a/P). This statement is at odds with the mechanism of polymer drag reduction and two other unique properties 
of polymer solutions, elastic turbulence and elasto-inertial turbulence. These three phenomena support de GenneiP's 
picture that drag reduction derives from two-way energy transfers between turbulent kinetic energy of the flow and 
elastic energy of polymers at small scales, resulting into the overall modification of the turbulence energy cascade. 
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The drag reducing mechanism is caused by an increase of the (extensional) viscosity in extensional upwash and 
downwash flows generated by quasi-streamwise vortices^-', thereby creating a negative torque on these near-wall 
vorticeiP. Dubief et alW demonstrated that polymers re-inject part of the energy accumulated in high-speed streaks, 
regions of locally high-speed flow in the near-wall region elongated in the direction of the how. In inertia-less flows 
with curved streamlines, Groisman & Steinberg^ demonstrated the existence of strong nonlinear mixing supported 
by elastic turbulence, a state of saturated dynamical interactions between stretched polymer molecules and the base 
flow that causes the stretching. 

The present paper focuses on the latest evidence of energy transfer from polymers to flow: the recently discovered 
elasto-inertial turbulence, hereafter referred to as ElT-S. EIT is a new state of small-scale turbulence driven by the 
interaction between elastic instabilities and the flow's inertia that has been observed over a wide range of Reynolds 
numbers, from subcritical to supercritical Reynolds numbers. EIT exists by either creating its own extensional flow 
patterns, as we will demonstrate here in subcritical channel flows, or by exploiting extensional flow topologies, as we 
shall observe in the wake of a hairpin vortex at higher Reynolds numbers. EIT offers an alternative description of 
MDR that has the merit of explaining the phenomenon of early turbulence^, which describes the onset of turbulence 
in the presence of diluted polymer additives at Reynolds numbers significantly smaller than in the absence of polymers. 
EIT also provides a unique insight on the possible coexistence of inverse energy cascade in three-dimensional flows. 



II. METHOD 



Channel flow simulations are performed in a cartesian domain, where x, y and z are the streamwise, wall-normal 
and spanwise directions, respectively. For a polymer solution, the flow transport equations are the conservation of 
mass, V • u = 0, where u is the velocity vector, and transport of momentum: 



d t u+ (u • V)u = -Vp 



Re 



1-P 
Re 



V T. 



(1) 



The Reynolds number is based on the bulk velocity Ut and the full channel height H = 2h, Re = U^H/v, The 
parameter j3 is the ratio of solvent viscosity to the zero-shear viscosity of the polymer solution and affects both the 
viscous stress and polymer stress terms in Eq. ([I]). The polymer stress tensor T is computed using the FENE-P 
(Finite Elastic Nonlinear Extensibility-Peterlin) models 
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where the tensor C is the local conformation tensor of the polymer solution and I is the unit tensor. The properties 
of the polymer solution are /3, the maximum polymer extension L, and the Weissenberg number Wi based on the 
solution relaxation time A and the flow time scale relevant to the dynamics of interest. Here Wi is based on the 
wall shear-rate 7 of the initial laminar flow at each Re, hence Wi = A7. The FENE-P model assumes that polymers 
may be represented by a pair of beads connected by a nonlinear spring defined by the end-to-end vector q. The 
conformation tensor is the phase-average of the tensorial product of the end-to-end vector q with itself, C = (q ® q) , 
whose transport equation is 



d f C + (u • V)C = C(V u) + (V u) T C 



(3) 



On the right-hand side of Eq. ([3]), the first two terms are responsible for the stretching of polymers by hydrodynamic 
forces, whereas the third term models the internal energy that tends to bring stretched polymers to their least energetic 
state (coiled). 

Eqs. (l|3) are solved using finite differences on a staggered grid and a semi-implicit time advancement scheme 
described elsewhere^. A series of simulations was carried out for Reynolds numbers ranging from 1000 to 6000. A 
thorough resolution study led us to choose a domain size of 10H x H x 5H with 256 x 151 x 256 computational nodes. 
All results discussed here have been verified on domains with a factor 2 in horizontal dimensions and resolution in 
each directions. The CFL number was set to 0.15 to guarantee the boundedness of C. 

The protocol for our simulations was designed to mimic the perturbed experimental setup of Samanta et a/P^ 
within the limitation inherent to the DNS boundary conditions. For any flow, Newtonian or polymeric, the initial 
flow and polymer fields are first equilibrated to the laminar state corresponding to the desired Reynolds number. A 
perturbation is then introduced over a short duration, in the form of blowing and suction velocity on both walls, over 
which white noise of prescribed intensity is introduced. The velocity pattern is periodic in x and z: 



v w (x,z,t) = H(t) 



A sin 
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FIG. 1. (a) Friction factor as a function of the Reynolds number for two Weissenberg numbers, Wi = 100 (•) and Wi = 700 (a), 

with L = 200 and /3 = 0.9. Lines indicate correlations for laminar ( , / = 12/ Re) and turbulent ( , / = 0.073i?e~ 1//4 ) 

Newtonian channel flow and for MDR ( , / = 0.42i?e -0 ' 55 ); Newtonian solutions are also included (□ ). (b) Mean velocity 

profiles normalized by viscous scales; polymer flows: o, Re — 1000, Wi = 100; □ , Re — 6000, Wi = 100; & , Re — 6000, 

Wi = 700; , Newtonian flow at Re = 6000. — , Virk's MDR velocity profile U + = 11.71n(y+) - 17. Inset: Difference 

between the analytical Newtonian laminar velocity profile and the velocity profile of polymer flow at Re = 1000, Wi = 100. 

(c) RMS of velocity fluctuations normalized by viscous scales: , streamwise velocity u ; , wall-normal velocity v'; 

, spanwise velocity w' for polymer flow at Re — 1000, Wi = 100. (d) Same as (c) for polymer flows at Re = 6000, 

Wi = 100 and Wi = 700, and for Newtonian flow at Re = 6000 (lines without symbols). For colors and symbols see (b). 



where A is the amplitude, L x and L z are the horizontal domain dimensions, and e{t) is the random noise. The total 
duration of the perturbation is 0.5h/Ub, of which the first and last 10% correspond to a gradual increase/decrease 
through a smooth step function T-L(t). Choosing A = 0.09C4 and the RMS of e at 0.005J7b causes the Newtonian flow 
to transition at Re = 6000. 

All polymeric simulations are performed with L = 200 and /3 — 0.9. Two Weissenberg numbers are considered, 
Wi = 100 and 700. The former is consistent with previous simulations of MDT^MHI The theory of Procaccia et aiP, 
based on infinite Re and Wi, motivates the second, as an exploration of the effects of very large elasticity, or Wi, on 
the flow. 



III. STATISTICAL DESCRIPTION OF THE POLYMERIC FLOWS 



A. Velocity statistics 

Figure [TJ summarizes the main statisticalproperties of our simulations, which were found to be in excellent agreement 
with experiments in a pipe flow^. Figure u\ a) shows the evolution of the friction factor as a function of the Reynolds 
number from Re — 1000 to 6000. The friction factor of polymeric flows departs from the laminar asymptote, with a 
slight drag increase, to transition smoothly to the MDR asymptote, where drag is reduced compared to turbulent flow. 
The effects of increasing elasticity is only noticeable beyond the junction between the laminar and MDR asymptotes 
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FIG. 2. (a) Mean polymer extension: o, Re = 1000, Wi = 100; □ , Re = 6000, Wi = 100; - Re = 6000, Wi = 700. (b) 
Polymer (same symbols as in a)) and Reynolds (lines) shear stresses normalized by the bulk velocity and channel half-height: 
, Re = 1000, Wi = 100; , Re = 6000, Wi = 100; - - , Re = 6000, Wi = 700; , Re = 6000, no polymers. 



of the friction factor at Re c = 1791. The lower Weissenberg number simulations reach a high drag reduced state 
(HDR 6 ), while the high Weissenberg number simulations achieve MDR. 

Figures[T|b-d) compare profiles of velocity statistics of the polymeric and Newtonian flows at two Reynolds numbers, 
Re — 1000 and 6000, representative of the two domains Re < Re c and Re > Re c . At low Reynolds numbers, there 
is no evidence of logarithmic behavior in the mean velocity profile at HDR or MDR, as discussed in a previous 
publication^. The mean velocity profile departs from the Newtonian Poiseuille solution for laminar flow (Figure [T]d) 
by only a few percents. In the second domain (Re > Re c ), mean velocity profiles differ significantly from Newtonian 
turbulence at the same Re. The highest Weissenberg simulation approaches Virk's MDR velocity profile^. When 
normalized by viscous flow scales Figure [ljc-d), the RMS of velocity fluctuations show a consistent picture at the two 
Reynolds and Weissenberg numbers considered here, that qualitatively resembles that of Newtonian wall-turbulence. 
The existence of a maximum in the RMS of streamwise velocity fluctuations suggest the presence of streaks in all 
flows. The intensity of these streaks and their thickness is obviously modulated by the elasticity of the solution, 
as obvious at Re = 6000 in Figure [ljd) from the shift of the maximum of RMS(u + ) from y + m 12 (Newtonian) 
to y + w 30 (Wi = 100) and y + ~ 50 (Wi = 700). The spatial distributions of spanwise velocity fluctuations for 
Re = 1000, Wi = 100 and Re — 6000 and Wi = 700 are qualitatively similar with a more significant decrease of 
the centerline RMS relative to its maximum than observed for the Newtonian and Re = 6000, Wi = 700 flows. In 
all polymer flows the distribution of the wall-normal velocity fluctuations is maximum at the centerline, which is a 
difference from Newtonian turbulence. 

The similarity between Re = 1000, Wi = 100 and Re — 6000, Wi — 700 suggests that the two flows create 
similar turbulent structures. We will later confirm that Re = 6000, Wi = 700 is indeed dominated by elasto-inertial 
turbulence. The difference, especially in RMS(w + ), with Re — 6000, Wi — 100 indicates the existence of a flow 
which is still under some influence of typical structures found in Newtonian wall turbulence. 



B. Polymer statistics 

a) shows the profiles of polymer stretch, denned as (tr(C)J /L, for the three simulations of interest. 
The mean polymer stretch is well below 50% of full extension at Wi = 100, indicating that the mechanism at play is 
not induced by coil-stretch transition, at least from the perspective of time-averaged transport equations. The most 
elastic solution has significantly larger mean stretch, yet its mean is well below 100%. In a purely laminar flow, i.e., a 
pure shear flow, the solution for the polymer stretch decreases monotonically and linearly due to the linear behavior 
of the shear. The emergence of a local maximum away from the wall at our lowest Reynolds number suggests the 
existence of extensional flow topologies, as such topologies are known to produce the largest polymer extensions in wall 
turbulence. Indeed, the local maximum of polymer at Re = 6000 arises from the interactions between polymers and 
turbulent structures in the buffer region, in particular vortices, which produce local extensional flows with dramatic 
effects on the polymer dynamics^. The possible existence of extensional flow topology at Re = 1000 is the motivation 
for our next section. 

Before studying the emergence of extensional flows, the Reynolds shear stress — uv and polymer stress (1 — (3)T xy / Re 
terms are compared (Figure [2b). We remind the reader that the two stresses, with the addition of the viscous stress 
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FIG. 3. Joint probability density functions of flow topology in the (Q a ,Ra) phase plot for polymeric flows at Re — 1000 

( ) and 6000 (-- - ) for Wi = 100 . Line describes D = 27/47^ + Q 3 a = 0, the vertical bounds of the four 

quadrants of flow topology. 



(not shown) , leads for the balance of stresses in a channel flow to: 



— 1 = —uv H 1 xv H 

da; V h) Re Re ay 



l ~Tl = ~ uv ^ 7^ T xy + 777 7777 ■ ( 5 ) 



For all simulations, the Reynolds shear stress is very small compared to the polymer shear stress. Supporting the 
proposition that flows at Re = 1000, Wi = 100 and Re = 6000, Wi — 700 are similar, the Reynolds shear stress is 
negligible in both cases, indicating the loss of a classical near-wall turbulence structure. The absence of Reynolds 
shear stress has noticeably been observed in the experiments of Warholic & HanrattjP^ for their MDR flow. Even 
at a slightly lower drag reduction, the polymer shear stress dominates the Reynolds shear stress, which is consistent 
with HDR ffowP. 

The statistical analysis of our polymer flow simulations shows turbulent states predominantly supported by polymer 
stresses. The similarity between the two cases Re — 1000, Wi = 100 and Re = 6000, Wi — 700 indicates that the 
state, EIT, observed at the lowest Reynolds numbers is likely to exist at higher Reynolds numbers, which suggests 
that EIT could be the asymptotic structural state of wall-bounded polymer flows. 



IV. TOPOLOGY OF EIT 

The flow topology is a critical component of the dynamics of polymers^, since it governs the stretching terms in the 
transport equation of the conformation tensor, Eq. ph, through the velocity gradient tensor Vu = djUi. We therefore 
apply the classical reduction of the flow into a joint p robab ility density function (jpdf) of the second Q a = —djiiidiUj/2 
and third R a = —djUidkUjdiUk/3 invariants of Vi i 16 * 17 l Based on the sign of the discriminant D — 27/AR^ + Q al 
quadrants I (R a > 0, D > 0) and II (R a < 0, D > 0) of Figure [3] define spiraling flows under compression and 
extension, respectively, and III (R a < 0, D < 0) and IV (R a > 0, D < 0) biaxial compressional and extensional flows, 
respectively. At Re = 6000, the jpdf contours exhibit an inverted teardrop shape common to many turbulent flows, 
in particular of Newtonian turbulent channel flow. The topology distribution for the lower Reynolds number flow is 
significantly different with a quasi-symmetry around both R a — and Q a = and the confirmation of the existence 
of biaxial extensional events (D < 0). 

Figure [4] shows the 3D structure of Q a with positive and negative isosurfaces of Q a for Re — 1000 and 6000. At the 
lower Reynolds number (Figure [4^,) , the structure of the Q a field is predominantly spanwise with trains of cylindrical 
structures of various scales. The smallest scale structures depicted by the threshold chosen in this figure are organized 
in trains of structures of alternating sign. Much larger scales are observed throughout the domain, with a typical 
spanwise dimensions comparable to the small scale structures. The typical streamwise dimension appears to be a 
couple of wavelength of the structures forming the trains of small scales. At Re = 6000 and moderate elasticity 
(Wi — 100), we chose the particular instant that is shown in Figure Qb) because the structure is characteristic of 
an active state^, or state of high drag (thus higher turbulence). No trains of alternating sign cylindrical structures 
of Q a may be observed. However, these trains reappear during hibernating stateiP^ (not shown for Wi — 100), as 
evident from Figure gc) for Re = 6000 and Wi = 700. Active states at Wi = 700 (Figure [§i) consist also of hairpin 
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FIG. 4. Instantaneous isosurfaces of the second invariant of the velocity gradient tensor for Q a — ±0.025 at Re — 1000, 
Wi = 100 (a); Q a = ±0.25 at Re = 6000, Wi = 100 (b); Q a = ±0.25 at Re = 6000, Wi = 700 during hibernating (c) and 
active (d) state. Red: positive Q a ; cyan: negative Q a . 



vortices, though much larger than for Wi — 100, and trains of small scale cylindrical structures of alternating sign of 
Q a . Some of these trains appear in the wake of heads of hairpin vortices and were also observed at higher Reynolds 
numberd^l Other train occurrences are also observed in regions of low turbulence, consistent with the hibernating 
state's picture. 

Local extensional flows are captured in contours of polymer stretch (Figure[5| in the form of thin sheets of locally high 
polymer stretch, tilted upwards and elongated in the flow direction. Superimposed to the contours of polymer stretch 
are the contours of positive and negative Q a , which reveal the existence of trains of circular regions of alternating 
Q a -sign associated to sheets of large polymer stretch. The intermittence of these sheets is clear for Wi = 100 at 
both Reynolds numbers, where the sheets of highly stretched polymers are surrounded by large regions of low stretch 
(Figures [5^,-b) . The smallest scales (cylindrical structures of positive and negative Q a ) cluster along the sheets of 
highly stretched polymers. The cross-section of these cylindrical structures are shaped like the outside contour of a 
dumbbell, or an ellipsoid squeezed along its shorter radius. The large scales of the Q a field appear also to alternate 
sign along the streamwise direction. These large scale structures seem to arise from the smallest scales maybe through 
a merger of the upper head (or most outward head) of the small scale dumbbell structures. This is only a speculation 
at this point which warrants further investigation but is not critical to the current study. 

In the active state of the Re = 6000, Wi = 100 simulation (Figure [5]d) , the trains of small scales are not present, 
at least in the cross-section considered here. At the highest Weissenberg number, highly stretched polymers occupy 
most of the domain, yet sheets are still present in the hibernating and active states, Figures [5^c) and (d), respectively. 
During the hibernating state, the structure of Q a exhibits the same dumbbell shape found at Re = 1000. The plane 
in Figure [5jcl) cuts through a hairpin vortex, illustrating one of the train of small-scale Q a structures spanning over 
— 1 < x/h <3 in the lower channel half. 

One important conclusion of this topological study is that EIT consists of cylindrical structures of Q a of alternating 
signs at least at two scales, small dumbbell-like cylinders and larger structures that seem to results from the merging 
of the dumbbell cylinders. Another important observation is that the features of EIT disappear when the flow is too 
turbulent or the polymer solution not elastic enough. EIT reappears during hibernation events or even during active 
events if the Weissenberg number is sufficiently large. This suggests that EIT is an asymptotic state that should occur 
when the elasticity of the solution can efficiently control and contain the growth of turbulence. 
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FIG. 5. Instantaneous contours of polymer extension in an x — y plane showing streamwise sheet-like regions. Re — 1000, 
Wi = 100 (a); Re = 6000, Wi = 100 (b) and Re = 6000, Wi = 700 for hibernating (c) and active (d) state. The continuous 
and dashed lines superimposed to the contours of polymer stretch represent isosurfaces of Q a : black, negative; white, positive. 
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FIG. 6. Profiles of the production of turbulent kinetic energy, Vh (lines), the negative contribution of dissipation rate of TKE, 
e (closed symbols) and transfer of energy between elastic energy and TKE, n e (open symbols). The sum of the volumetric 

integral of these term results in the energy balance found in Eq. Q. Re = 1000, Wi = 100: , o , • ; Re = 6000, 

Wi = 100: , □ , ■ : Re = 6000, Wi = 700: , a , . 



V. SIMILARITY BETWEEN EIT AND 2D TURBULENCE 

A. Energy transfer between polymers and turbulence 

The only evidence presented so far for the sustainability of EIT through energy transfer from polymers to turbulence 
is the very existence of EIT at a Reynolds number much lower than the critical Reynolds number for Newtonian flows. 
Following Dallas et aZp3 we use the integral of the budget of turbulent kinetic energy (TKE): 

V k dV- [ edV- [ U e dV = 0, (6) 
v Jv Jv 

where V is the volume of the channel. The three terms involved are the production of TKE, the dissipation rate of 
TKE and the energy transfer between polymers and TKE, which are respectively defined as: 

v h = -«v? ( 7a ) 

ay 

e = ^&7W (7b) 
Re 



n e = —^S' : T' (7c) 
Re 

The tensor S' is the fluctuating strain rate tensor S' = (Vu' + Vu' T )/2 and T' is the fluctuating polymer stress 
tensor. Figure [6] shows the profile across one half of the channel of the three terms of interest to for polymer flows 
at Re = 1000 and 6000. The wall-normal coordinate is normalized by viscous scales. In Eq. Q, the production 
and dissipation rate of TKE terms are always positive. The energy transfer between polymer and turbulence divides 
the flow in two regions. In a thin region attached to the wall (y + < 7, polymers supplements the dissipation rate 
of TKE. Beyond this "polymeric viscous sublayer" , polymers support and even create, for Re = 1000, turbulence 
with transfers from elastic to turbulent kinetic energy. Lending further credence to our proposition that EIT is the 
asymptotic state of polymer drag reduction, the production of TKE in the highest Weissenberg and Reynolds number 
flow is overwhelmingly dominated by — II e , in a similar fashion as for the Re — 1000 flow. 

Figure [6] brings further support to de Gennes's theorj0 of exchange of energy between elastic energy and TKE. 
A similar point was made by Dallas et alW^ using DNS as well. The picture of EIT that emerges so far is one of 
extraction of TKE through the damping mechanism of coherent structures^ and injection of elastic energy into TKE 
through elastic instabilities caused by the formation of thin sheets of polymer stretch, as highlighted in the previous 
section and Figure [5| 
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FIG. 7. (a) Spectra of turbulent kinetic energy at y + — 15. (b) Spectra of Q a at the same location. : Re = 1000, 

Wi = 100; Re = 6000, Wi = 100; Re = 6000, Wi = 700. Inset: Spectrum of Q a for Newtonian turbulence at 

Re — 6000, same location. 
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FIG. 8. Spatial correlations in the spanwise direction of streamwise (a) and wall-normal (b) velocity fluctuations at y + = 15. 
Line captions for polymeric flow is the same as in Figure [7] , Re = 6000 Newtonian; 



B. The 2D turbulence analogy 

Equipped with the knowledge of energy transfer from polymers to TKE, we now turn to spectral analysis for a 
discussion on the mechanism of such transfer. An interesting characteristic of the spectral distribution of TKE in 
EIT flow is the observation of k a behavior at high wavenumbers. Figure [7][a) shows longitudinal ID spectra of TKE 
at y + — 15 for all polymeric simulations. The location was chosen for its good approximation of the maximum 
energy transfer of elastic energy to TKE for all simulations. In spite of the low Reynolds numbers considered here, 
a range of wavenumbers, between k x w 4 to 10 for Re = 1000 and k x « 8 to 50 for Re = 6000, show a power 
decay that fits well /c~ 14 / 3 . We checked the robustness of the power law behavior by plotting the exponent derived 
from a(k x ) = {k x / E(k x )){dE / dk x ) and confirming the existence of a plateau of a(k x ) in the range of wavenumbers of 
interest (not shown). The —14/3 exponent bears, at this point, no theoretical significance; the best fit was found to be 
around 4.6. In the two EIT flows, the higher wavenumbers depart from —14/3 and seem closer to —11/3 (not shown). 
Nonetheless the important observation here is that the power law decay is between —3 and —5. Figure[7](b) correlates 
the lower bound wavenumber of the —14/3 region with the maximum spectral contribution to the fluctuations of Q a . 
This energetic peak is dramatically more marked in polymeric flows than in Newtonian turbulence, and corresponds 
to the trains of cylindrical structures observed in the instantaneous visualizations shown in Figures [4] and [5] For 
reference, fc -5 / 3 is also fitted to the lowest wavenumbers. The reasonable visible agreement is more an effect of the 
vertical scale of the plot rather than an actual —5/3 behavior, yet this reference will help the present discussion. 

The spectral behavior described here provides circumstantial evidence of a 2D-likc dynamics of polymeric flows at 
or close to EIT. Indeed, 2D flows are known for their clearly observed, documented and even theoretically described 
inverse energy cascade^. In fact, the best analogy between our simulations and 2D turbulence is found in numerical 



10 



experiments where friction is introduced^. In such experiments, the flow is governed by the following incompressible, 
2D Navier Stokes equation: 

d t u+ (u • V)u = -Vp + ;A7 2 u - v f u + f , (8) 

where i>f is the friction coefficient and f is a discrete random forcing at a given intermediate wavenumber kd- These 
flows are solved in periodic domains. Such flows experiences an inverse energy cascade as welP^, where energy cascades 
toward the largest scales (k < kd) with an energy spectrum of the form: 

E(k) = C K e$k~% for k < k d . (9) 

At larger wavenumbers, energy is dissipated through an enstrophy cascade, where the enstrophy is the integral of 
the square of the vorticity. The influence of the friction coefficient is felt in the exponent of the spectral decay of 
enstrophy cascade which corresponds to 

E(k) oc fc" 3 ~« , (10) 

where < £ < 2, in agreement for the slopes found in our polymeric simulations. The analogy is further supported by 
the significant increase of the correlation length in the spanwise direction as shown in Figure [8] for the streamwise and 
wall-normal velocity fluctuations. For the lowest Reynolds number simulation and the highest elasticity at Re = 6000, 
the correlation length is dramatically increased. The flow does not experience however a reduction in dimensionality 
as dramatic as the ones observed in strongly rotational flowPl for instance. It is however striking to note that, 
even at the lowest Reynolds number, the structure depicted by the correlation of the streamwise velocity fluctuations 
indicate the existence of streak-like structures, yet much wider than for the Newtonian flow. This was confirmed by 
instantaneous visualizations of wall shear stress (not shown). 

The concept of inverse energy cascade does however makes sense in the context of our polymeric flows. There is 
an obvious resemblance between the polymeric flow momentum equation Eq. ([!]) and the 2D equation Eq. ^ in the 
dynamics at play, where the divergence of polymer stress embodies both the friction term and the injection of energy 
term. Since the same term provides large scale dissipation (drag reduction mechanism) and energy injection, one may 
expect differences between the two dynamical systems, which could explain the coexistence of the -14/3 and -11/3 in 
EIT flows. 

The elastic instabilities captured through isosurfaces and contours of Q a (Figures |1J [5] and [7)3) occur at intermediate 
wavenumbers close to the lower bound of the -14/3 range shown in Figure [7j Whether such instabilities excite or 
directly create larger scales is a question we cannot answer with the present analysis. One possible approach that 
could provide valuable insights in the future is the method proposed by Biferale et al^Q which demonstrated the 
existence of inverse energy cascade in 3D isotropic turbulence. Incidentally, Biferale et al. demonstrated that "2D 
and 3D properties naturally coexist in all flows in nature" . The interesting conjecture derived from their study and 
our observations here is that a flow that is locally two-dimensional may exhibit a behavior akin to that observed in 
true 2D turbulence, in particular an inverse energy cascade. 



C. An instantaneous picture of transfer between elastic and turbulent kinetic energies 

Assuming that the proposed analogy between EIT and 2D turbulence holds, the energy transfer term between 
polymers and TKE can be interpreted as follows. TKE is a large scale quantity, whereas the organization of polymers 
in sheets strongly suggests that polymers are intermediate or small scales, as shown by the Q a structures they create. 
Indeed the thickness of the sheets of highly stretched polymers appear to govern EIT, as evident from the clustering 
of Q a structures on these sheets. Defining the volumetric integral as (•) = J v »dV, the rate of change of TKE is thus 

= {Vk) (e) _ ( n e) (ii) 

and — (e) is always negative thus a sink of TKE. The scalar product between tensors T' and S' is effectively the 
rate of work done by the strain rate on the polymer stress. When fluctuations of polymer and strain tensors align, 
T' : S' > 0, the work is a sink of turbulent kinetic energy, consistent with the increase of polymer stress in sheets 
(extensional flows) which triggers large increases in extensional viscosity. Following our analogy, this process would 
describe a forward energy cascade. When polymer stress fluctuations oppose flow deformation fluctuations, — IL e 
becomes positive and contributes directly to the TKE, which amounts to a backward energy cascade. The interplay 
between — II e and coherent structures isolated by isosurfaces of Qa is shown in Figure[9j The region of forward cascades 
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FIG. 9. Spatial correlations of streamwise (a) and wall-normal (b) velocity fluctuations in the spanwise direction. 
Re = 6000 Newtonian; 



These regions also 
9|j for < x/h < 4, 



are sheets, that one can visually correlate to sheets of highly-stretched polymers from Figure p 
appear to be the cause of the narrow part of the dumbbell shape of the Q a -isosurfaces (Figure 
for instance). The backward energy cascade surrounds the sheets of forward cascade and contain the largest scales of 
Q a . Note that in the active state depicted in Figure [9jci) , one can observe the wake of the head of a hairpin vortex 
where the flow is extensional and triggers EIT. The distribution of backward and forward energy cascade is evident 
yet very complex and warrants future studies. Even at the lowest drag reduction studied here, there is still significant 
amount of backward energy cascade. 
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FIG. 10. Conceptual skecth of the mechanism of EIT. 



VI. PROPOSED MECHANISM OF EIT 



The dynamic of EIT was found to be consistent across the two regimes shown in Figure [TJa) (Re < Re c and 
Re > Re c ), and can be best described by (i) taking the divergence of Eq. 0, which yields the Poisson equation for 
pressure in a viscoclastic flow: 



V 2 p = 2Q a - V^ V ' ( V ' T )' ( 12 ) 
Re 

and (ii) considering the hyperbolic nature of the transport equation of the conformation tensor, Eq. caused by 
the absence of diffusion^. For Re < Re Cl a perfectly laminar flow stretches polymers through the action of shear 
(stretching term in Eq. [3]). The introduction of small perturbations into the flow excites the unstable nature of the 
nonlinear advection term [u- V)C, resulting in the formation of sheets or cliffs of polymer stretch akin to cliffs of scalar 
concentration observed in the turbulent transport of low-diffusivity passive scalar^. This behavior is obvious when 
comparing Figure [3](b-c) with LIF (Laser Induced Fluorescein) images of fluorescein dye concentration in polymer 
drag-reduced wall bounded fiows^J. The sheets of high polymer stretch, hosting a significant increase in extensional 
viscosity, create a strong local anisotropy, with a formation of local low-speed jet-like flow. The response of the flow 



is through pressure (Eq. 12 1, whose role is to redistribute energy across components of momentum, resulting in the 
formation of waves, or trains of alternating rotational and straining motions as shown by the Q a isosurfaces. The 
mechanism shares some similarity with the Kelvin-Helmholtz instability, except that the thickness of these sheets is 
too close to the Kolmogorov scale (smallest scale of turbulence) for vortices to be created. Once triggered, EIT is 
self-sustained since the elastic instability creates the very velocity fluctuations it feeds upon. Interestingly, Samanta 
et al^H and Dubief et aZP^ show that the phenomenon of EIT is not confined to low Reynolds numbers, unlike 
elastic turbulence^. We also presented circumstantial evidences that EIT shares some interesting similarities with 2D 
turbulence, which could provide a reasonable framework for further studies. The mechanism of EIT is summarized in 
Eq. JlOl). 



VII. CONCLUSION 



Elasto-inertial turbulence offers a new perspective on polymer drag reduction. First, it provides support to 
DeGennes^ theory of energy transfers between polymers and flow. Second, EIT allows us to consider the possi- 
ble structure of MDR for very large elasticity (Wi — > oo): The sheer magnitude of extensional viscosity is likely to 
prevent the emergence of any vortical structures, thus leaving MDR to be sustained by near-wall spanwise structures 
similar to the ones observed at low Reynolds numbers (Figure [3ji) . The absence of vortices, however, does not signify 
the absence of streaks. Indeed, there is no apparent means in our current understanding of EIT to damp streaks. Our 
study suggests that, given the right amount of injected energy, the streak instability can be triggered and sustained 
at subcritical Reynolds numbers. 

The asymptotic state of polymer drag reduction in the absence, or near absence, of vortices should therefore be 
driven by the nonlinear transport of polymer stretch, which resembles the transport of a high Schmidt number passive 
scalar, and the response of the flow to a sheet-like, strongly anisotropic field of effective viscosity governed by the 
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extensional viscosity of polymers. As discussed by Dubief el al\ 19 ^ 27 \ the flow is therefore stuck in a transitional state, 
specifically the stage of breakdown of nonlinear flow instabilities, which does not support a logarithmic mean velocity 
profile 28 . Unfortunately, a major obstacle in the derivation of a low Reynolds (or possibly high Reynolds) number 
correction of Virk's log law is the lack of theoretical understanding of high Schmidt number active scalar transport 
in nonlinear anisotropic flows. Our analogy with 2D turbulence speculates that the mechanism by which turbulence 
is supported may be that of an inverse energy cascade. The potential broader impact of this discussion could be that 
inverse energy cascade in 3D turbulent flows may be related to the local reduction dimensionality of the flow. Another 
perspective is also the possibility to develop a theoretical approach using a 2D or quasi 2D flows. 

Future research should therefore focus on the dynamics of active scalars and the backscatter of energy caused by 
polymer dynamics. Another interesting future research direction is whether the action of the polymers in the flow 
is catalytic or one of direct energy exchange between polymer and flow instabilities. At low Reynolds numbers, the 
process appears to be a direct energy exchange; however, certain aspects of EIT at high Reynolds numbers may be 
catalytic, e.g. the trigger of elastic instability through bypass transition. 
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